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1. Introduction 



Curvature in relativistic spacetimes corresponds to tidal forces in Newtonian 
mechanics, but curvature effects yield more precise information about physical 
phenomena. In particular, general relativity should provide corrections to cal- 
culations based on Newtonian physics for the statistical mechanical behavior 
of a gas subject to a gravitational field. Such an investigation was undertaken 
in pQ, where an approximate Helmholtz Free Energy function was derived, us- 
ing statistical mechanical methods, for an ideal gas in orbit around the central 
mass in Schwarzschild spacetime. For a partial listing and summary of research 
in relativistic statistical mechanics, we refer the reader to the introduction in pQ. 

In this paper, we consider an ideal gas enclosed in a container with coordinates 
at rest with respect to an observer whose four-velocity is a timelike Killing vec- 
tor, and we assume further that the gas particles, along with the container, do 
not signicantly alter the background metric. It should be noted that the term 
"ideal gas" is somewhat misleading in the context of general relativity. This is 
because a volume of gas subject to no forces is still affected by the curvature of 
space-time, and this corresponds to a Newtonian gas subject to gravitational, 
tidal, and in some instances "centrifugal forces," but otherwise "ideal." We 
provide details for this association in Sections 6 and 8. 

The significance of the frame of reference in relativistic statistical mechanics 
was discussed in Section 5 of [5] . There, in the context of Minkowski spacetime, 
it was argued that there is no Lorentz invariant thermal state for an ideal gas 
system. A thermal state is in equilibrium only in a preferred Lorentz frame, 
and therefore breaks Lorentz invariance. In the Lorentz frame in which the 
container of the gas is at rest, the Lorentz covariant canonical distribution den- 
sity for an ideal gas reduces to a product of terms of the form expftpo, where 
Po is the negative of the energy of a particle (in that frame) and f3 is inverse 
temperature. In that frame, the relativistic statistical mechanical behavior may 
be easily compared to its non relativistic analog. 

In the context of general relativity our requirement that the spacetime possess 
a timelike Killing field is physically natural. The four-velocity Killing vector for 
the container determines a time coordinate for the system and thus an energy 
component of the four-momentum. As a consequence, comparisons to corre- 
sponding non relativistic statistical mechanical formulas are possible. Moreover, 
this choice of four-velocity forces the gas system to evolve in such a way that 
the geometry of spacetime is unchanging along its worldsurface (since the Lie 
derivative of the metric vanishes in that direction) , and thus it is plausible on 
physical grounds that the particle system will reach equilibrium. 

In order to take full advantage of these features, we introduce, in Section 2, a 
non rotating, orthonormal "Fermi- Walker-Killing" coordinate system in which 
the gas container is at rest, and for which the Killing vector serves as the time 
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axis. In Section 3, we identify the appropriate volume form on phase space, and 
then prove Liouville's Theorem in Section 4, improving the local, approximate, 
version of Liouville's theorem established in 1 . In Section 5, we derive an exact 
relativistic formula for the Helmholtz free energy of an ideal gas. In contrast to 
PP, the Boltzmann factor in our partition function is an invariant scalar. In Sec- 
tion 6 we systematically find Newtonian limits for spacetimes and Killing fields 
satisfying general conditions. We compare the exact relativistic expressions to 
their Newtonian analogs, derived both independently and as Newtonian limits 
of our formula in Section 8. We also find in Section 7 the relativistic thermody- 
namic equation of state. Concluding remarks are given in Section 9. 

2. Fermi- Walker-Killing Coordinates 

Throughout we use the sign conventions of Misner, Thorne and Wheeler [3]. A 
vector field V in a spacetime M. is said to be Fermi- Walker transported along 
a timelike path a if V satisfies the Fermi- Walker equations given by, 

F4D = v 5 y a + n a / = o. (i) 

Here u is the four- velocity along a (i.e., the unit tangent vector), f2 a o = a a ug — 
u a ag, and a a is the four-acceleration. As usual Greek indices run over 0, 1, 2, 3 
and lower case latin over 1,2,3. It is well-known and easily verified that 
Fu{u) = 0, and if vector fields V and W are Fermi- Walker transported along 
a, the scalar product V^Wg is invariant along a. Thus, a tetrad of vectors, 
Fermi- Walker transported along a and orthonormal at one point on a, is nec- 
essarily orthonormal at all points on the path. Moreover, such tetrads may be 
constructed so that one of the orthonormal vectors is the tangent vector u. 

Let <t(t) denote the parameterization of a by proper time r, and let eo(r), 
ci(t), 62(7"), e3(r) be an orthonormal Fermi- Walker transported tetrad along tr, 
with e = u. The Fermi- Walker coordinates 3 relative to this tetrad 

on a are given by, 

x ° ( ex Pa(r)( AJe .?( r )) = T 

(2) 

x K (exp CT(T) (A^(T)) =X K , 

where exponential map, exp p ({;), denotes the evaluation at affine parameter 1 
of the geodesic starting at the point p in M, with initial derivative v, and it 
is assumed that the A J are sufficiently small so that the exponential maps in 
Eq.([2|) are defined. From the theory of differential equations, a solution to the 
geodesic equations depends smoothly on its initial data so it follows from Eq. ([2]) 
that Fermi- Walker coordinates are smooth. Moreover, it follows from [4] that 
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there exists a neighborhood of a on which the map X = (x°, x 1 , x 2 , x 3 ) is a 
diffeomorphism onto an open set in R 4 and hence a coordinate chart. 

Assume now that J? is a timelike Killing vector field in a neighborhood of a 
and the tangent vector to a is K, i.e., eo = u = K. The vector field K 
is the infinitesmal generator of a local one-parameter group <f> s of diffeomor- 
phisms. The function <j) s is the flow with tangent vector K at each point, and 
<t(t) = 4> T {0, 0, 0, 0) (where (0, 0, 0, 0) is the origin in Fermi- Walker coordinates). 

Define a diffeomorphism X^ 1 from a sufficiently small neighborhood of the 
origin in 1R 4 to a neighborhood U of cr(0) in M. by 

X~\x*, = ^opf-HM 1 ,^* 3 )) (3) 

Then X = < 3 ) is a coordinate system on U which we refer to as 

Fermi- Walker-Killing coordinates. The following properties are readily verified: 



_d_ 

dx, 1 



- A 
— = K 

dx° 



x°=0 



(4) 



Thus, in Fermi- Walker-Killing coordinates, the time coordinate x° is the pa- 
rameter of the flow generated by the Killing field K with initial positions of the 
form (0 3 ) in Fermi- Walker coordinates. 



In the sequel, we will make use of both Fermi- Walker and Fermi- Walker-Killing 
coordinates. In the case of Fermi- Walker coordinates we designate momentum 
form coordinates as p a and in the case of Fermi- Walker-Killing coordinates, the 
momentum form coordinates will be designated as p a . Thus, coordinates of 
the cotangent bundle of M will be represented as {x a ,pp} or {x a ,pp}. The 
metric components in Fermi- Walker coordinates are designated as g a /3 and in 
Fermi- Walker-Killing coordinates as g a p- Note that g a p does not depend on x Q 
because of Eqs.([4|). 

3. Phase Space 

Consider an observer following a timelike path a with four-velocity K which 
is a Killing vector in Ai, and assume that there is a gas of non interacting 
test particles enclosed in a container whose coordinates are at rest with respect 
to this observer. We assume further that neither the gas nor the container 
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significantly contribute to the gravitational field, so that the metric on M is 
effectively unaltered by their presence. In a neighborhood of a, spacctime points 
may be charted using Fermi- Walker coordinates or Fermi- Walker-Killing coor- 
dinates (described in the previous section), whose origin is <x(0). 

The Hamiltonian for the system of N ideal gas particles is given by 

N I 
H = 2jflj, where H, = -gf p p Ia Pip, (5) 

1=1 

and where Hj is the Hamiltonian for the Jth particle. With the Hamiltonian 
expressed in terms of Fermi- Walker-Killing coordinates, the natural time coor- 
dinate, and parameter for the dynamics, of the N particle system is 5°. We 
consider the dynamics in the following section. In the present section, we derive 
the basic 6-form for the relevant phase space. 

For simplicity we omit the subscript / when it is clear that we are referring to 
a one-particle Hamiltonian. The state of a single particle consists of its four 
spacetime coordinates together with its four-momentum coordinates {x a ,pp}, 
but the one-particle Hamiltonian satisfies, 

H: = \g aP P a p P = -^f, (6) 

where, for timelikc geodesies, m is the proper mass of the particle. Eq.([(j]) may 
be rewritten as, 



_ -g°% + y/(g°3g 0k - g ot) g 3k )P 3 Pk - f°m 2 c 2 

Po = ^ • (7) 

This allows us to reduce the dimension of the one particle state space to seven. 
Let P be the sub bundle of the cotangent bundle of spacetime with each three 
dimensional fiber determined by (O . The proof of the following proposition is 
given in pQ . 



Proposition 1 The volume 7-form Cj on P given by, 

Co = -^-dx A dx 1 A dx 2 A dx 3 A dpi A dp 2 A dp 3 , (8) 
is invariant under all coordinate transformations. 

Phase space for a single particle is determined by the space slice, at fixed time 
coordinate, orthogonal to the observer's four-velocity (along cr(r)), along with 
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the associated momentum coordinates, i.e., the cotangent bundle of the space- 
slice at fixed time coordinate x . The appropriate volume form is given by the 
interior product i(md / dT)Co of the four-momentum vector md/dr = p a d/dx a 
with lu. Then, 

mi{d/dT)uj = dx 1 A dx 2 A dx 3 A dp 1 A dp 2 A dp 3 + (dx° A VO, (9) 

where ip is a five-form. Since dx° = on vectors on phase space (with fixed 
time coordinate the restriction of mi(d/dT)uj to the one-particle (six- 
dimensional) phase space is, 

mi(d/dT)Cj = dx 1 A dx 2 A dx 3 A dpi A dp 2 A dp 3 . (10) 

It follows from Proposition[T]that the 6-form given by Eq. lfTO]) is invariant under 
coordinate changes of the space variables with x° fixed. Physically this means 
that the calculations that follow below are independent of the choice of the or- 
thonormal triad ei(r), e2(r), & 3 {t). 

4. A Liouville Theorem 

The phase space volume form for particles is the product of N copies of 
Eq. fTU)) , one copy for each particle. We show that this measure satisfies a Li- 
ouville theorem. The phase space, defined in terms of Fermi- Walker-Killing 
coordinates, is 6./V dimensional, but for ease of notation, consider instead an 
ideal gas consisting of N particles in a one-dimensional box; the generaliza- 
tion to additional degrees of freedom is straightforward. Label the coordinate, 
momentum, and proper time of the Jth particle, (I = 1, . . . , N) by x^pi^r^ 
respectively. We assume that t, — Tj(x a ), where x° is the time coordinate in the 
Fermi- Walker-Killing coordinate system. Our argument requires only that the 
function ^(a; ) be a smooth function of x° (in the interval between collisions). 

As noted previously, the Hamiltonian, H,, governing the dynamics of a sin- 
gle particle given by Eq. © is independent of x°. We may thus define the 
Hamiltonian vector field, vh, as follows, 



G 



( Si \ 

Pi 

xn 
\ Pn / 



( —L \ 

dx° 

dpi 
dx° 



dx 



N 



V 



dx° 

dpN 
dx° 



I dxi dr\ \ 
c?Ti (£r° 

dpi dr\ 
dri dx° 



dx 7v rfrjv 
dpN dr^ 



V 



(11) 



drjv dx° 

Liouville's theorem, which states that the Hamiltonian system is measure-preserving, 
is established from the following calculation: 



,. „ dxi dpi 

div v H = ■—- + -—- + 

OXi dpi 

d ( dxi 
dxi \dn 
d 2 H 



dx 



N 



dp 



N 



+ 



dXM 

d ( dpi 
dpi 



d 2 H 



dxidpi dpidxi 



dri 
dri 
~dlP 



dpN 
dri 
dx° 



+ ••• + 



dr N 
dx ' 



dr N 
dx° 



0. 



(12) 



It now follows that the iV-fold product of 6-forms of the form, dx 1 A dx 2 A dx 3 A 
dpi A dp2 A dp3, is invariant under the Hamiltonian flow. 

5. The Canonical Free Energy 

Non relativistic statistical mechanics expresses time averages of observables as 
phase space averages (in which time is no longer a parameter). In this sec- 
tion, we define and calculate the canonical partition function for a relativistic 
particle system at rest with respect to the Fermi- Walker-Killing coordinate sys- 
tem, developed in previous two sections. For this purpose, we will identify the 
phase space of a single particle as the cotangent space of the space-slice at a 
fixed arbitrary time 5°, and for convenience we take x° = 0. This is the col- 
lection of simultaneous events at time a; = together with associated momenta. 

The four velocity of the observer following the path ct(t) is K. The components 
K a of K in Fermi- Walker coordinates, and K a in Fermi- Walker-Killing coordi- 
nates, on (t(t) are (1,0,0,0). The observer, ct(t), then measures the energy of 
a particle with four- momentum p to be —K a p a = —p in Fermi- Walker-Killing 
coordinates, or —K a p a = — p in Fermi- Walker coordinates. Since the scalar 
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product is invariant under changes of coordinates, and K is a Killing vector 
field, —K a p a = —K a p a is constant along any geodesic tangent to the momen- 
tum vector in the container of gas. Thus, it is natural to define the energy 
of a particle with mass m and momentum p a to be —K a p a , but we modify 
this definition of energy by the addition of the constant —mc 2 , for the purpose 
of later comparison of relativistic statistical mechanics and classical statistical 
mechanics (which traditionally does not take account of the rest mass energy). 
The energy of a particle is then —K a p a — mc 2 , and the total energy of the 
particles in the container is the sum of the energies over all N particles. 

For what follows, note that by virtue of Eq.©, the energy —K a p a — mc 2 of 
a particle is uniquely determined by K and the space components pi of the 
momentum form. 

Let, 

dx = dx dx dx , dp — dp\dp2dpz, 

and 

d5i N dp N = dx\dx\dx\dpiidpi2dpiz . . . dx]^dx 2 v dx%dpNidpN2dpN3- (13) 

By Liouville's theorem, established in the previous section, and because particle 
energies, —K a p a , are constants of the motions along geodesies, it follows that 
for any fixed /3, 

e PT.^K«p a+ mc-) ld -N df> N (w) 

is an invariant measure on phase space for the N particle system (the index / 
in Eq. (fTi|) labels the particles in the container, as in Eq.©). In fact it is not 
difficult to verify that, 

£d/dx» (e^ Rapa+mc ^d5tdp\ = C K (e^ R "P" +mc ^dxdpj = 0, (15) 

where C denotes Lie derivative. Analogous to the development of nonrelativistic 
statistical mechanics, this invariance is a partial justification for the choice of 
HU) as phase space measure. We assume that our particle system is in contact 
with a heat bath, and consequently is in thermal equilibrium. Following the 
usual convention, let (3 = 1/fcT, where T is the temperature of the gas in some 
volume V, and k is Boltzmann's constant. 

Remark 5.1 An alternative convention for temperature may be used here, con- 
sistent with one discussed by Tolman in the context of thermodynamics [S] (see 
also [6]). At any position in the gas container, one may identify an observer with 
four- velocity fC/||A"||. The energy of a particle with momentum p measured by 
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that observer is — K a p a /\\ K\\ . Then, defining a position dependent equilibrium 
temperature T by T = r/||if|| = T '/ 'y '—.goo coincides with Tolman's formula 
and results in the same expressions we find below after normalizing by f3mc 2 . 

Now, for any fixed time coordinate, x°, we may define the canonical partition 
function for the gas in terms of Eq. (fT4|) . For convenience of calculation, we 
choose x° = r = 0. It then follows from Eq.Q that the phase space measure at 
the zero time coordinate is equal to the analogous expression in the (unbarred) 
Fermi- Walker coordinates. 

Thus, we define the canonical partition function for the ideal gas by the following 
integral expression in Fermi- Walker coordinates: 

N\(2ttH) 3N J v n J R 3n v k 1 

with notation defined analogously to (jT3|) . 



The volume integrals in Eq. (|16|) have limits of integration determined by the 
range of Fermi- Walker coordinates, x 1 , x 2 , x 3 , that define the volume of the gas, 
and the momentum integrals may be calculated explicitly. To carry out that 
calculation, consider the single-particle canonical partition function, 



1 



(2nh) 



>V 



and define, 



e 0K« Pa 



dp. (18) 



Since K and p are both timelike, 



^ Q p Q = -||if||||p||cosh X , (19) 

where x is the hyperbolic angle between the two vectors. In Eq. (|19p . ||p|| = 
mc, and we may write, \\K\\ — ac, where a = a(x 1 , x 2 , x 3 ) is a dimensionless 
function of the space coordinates. To evaluate Eq. fTS)) , we use the coordinate 
transformation on R 3 given by, 

Pi = mcsinhxcosfl, 

P2 = mc sinh \ sin 9 sin <fi, (20) 
ps = mc sinh \ sin 9 cos <j>, 

where ^ cj) < 2ir, ^ 9 < it, and x ^ 0. The volume element is 
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dpidp2dp3 — (mcf sinh 2 \ cosh \ sin 9 d\ dO d(f>. (21) 
Eq. (fT5|) then becomes, 

= 47r(mc 



o Jo 



/•oo 

) 3 / e - Q/3mc2coshx sinh 2 x cosh X rfx (22) 
Jo 

4 r ^ f° - 7a /~2 T W 47r(mc) 3 Jf 2 ( 7 ) 
47r(mc) / e w y\/y — 1 ay = , 



where y = coshx,7 = j(x , x 2 , a; 3 ) = a(ar, x 2 , i 3 )/3mc 2 , and -^2(7) is the 
modified Bessel function of the second kind. Substituting the above result for 
z p back in the expression for the single-particle canonical partition function, 
Ea.([T7|). we obtain 

_ 4tt(toc) 3 gme 2 /■ jr 2 ( 7 ) 

(2^) 3 iy 7 

Thus, Eq. (fH)]) may be rewritten as, 



-dx. (23) 



• r . q ts n \ 



(24) 



V2tt«^ J v 7 

Following the canonical statistical mechanical prescription, the connection to 
thermodynamics is given by, 



F(f3,V,N) = -~lnZ, (25) 

where F(f3, V, N) is the Helmholtz free energy of the gas. Note that the free 
energy given by Eq. (|25p is a function not merely of the volume of the container, 
but its actual shape via Eq. (|24p . This is an unavoidable consequence of the non 
uniformity of the gravitational field. 

In the special case of an observer following a timclike geodesic in Minkowski 
spacetime, the integrand in Eq. (l24|l is constant, and the free energy may be 
evaluated explicitly. The result is well-known, see e.g., [7]. In that case the 
observer has a four-velocity in Fermi coordinates given by K — (1, 0, 0, 0). With 
these components the vector field K is a Killing field, and it is clear that a(A) = 



X^-K a K a = 1. Thus, from Eq.(j23]), 

_ 0mc 2 47r(mc) 3 K 2 ((3mc 2 



(27rft) 3 (3mc 2 V: ^ 



and, 
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F03.V.JV) = - Win ^ - (27rfi)3 ^ yj . (27) 

consistent with [7] . The presence of gravity generalizes the special relativistic 
result so that 7(x 1 , x 2 , x 3 ) = a(x 1 ,x 2 ,x 3 )f3mc 2 becomes a function of spatial 
coordinates, rather than a constant. 

6. Newtonian Limit of the Relativistic Free Energy 

For a timelike path a tangent to a Killing vector if in a given spacetime, we 
have defined via Eqs. (f24|) and (j25|) the relativistic Hclmholtz free energy of 
a container of an ideal gas determined by a. With appropriate assumptions 
(described below), associated to such a statistical mechanical system is an anal- 
ogous Newtonian (i.e., nonrelativistic) expression for the Helmholtz free energy 
of the gas. This associated Newtonian Helmholtz free energy is the Newtonian 
limit of Eq. (f2"5|) . i.e., the limit as c — > oo, where c is the speed of light. 

The asymptotic behavior of ^2(7) for large argument is given by (see [5]), 

^2(7) ~ J^- e" 7 - (28) 
V 2 7 

and therefore, 



g/3mc (1 — a) 



2ti7i 2 /V Jv a 3 / 2 
In order calculate the Newtonian limit, we let 



(29) 



X = - c , (30) 

and we assume that a = a(\) is a smooth function of A (here and below we 
suppress the dependence of a on x , x 2 , x 3 ). Under fairly general conditions 
(see the following section), a(0) = 1 and a'(0) = 0, where the prime denotes 
differentiation with respect to A. In that case we may write, 

a(X) = l + ia"(0)A 2 + O(A 3 ) (31) 



Combining Eqs. ([29j) and (|3lj) . letting A — ► , and denoting the dependence of 
the one particle partition function on A by z\, yields, 
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/ m \I f e -i/3m(a"(0)+O(A)) 

A— +0+ ZA = (, 2^/3 ) A !™ + / y (i + I a "( )A2 + O(A3))3/2 rfx 

3 ^ ) 

e -i/W' ( o) rfx = ZNcwL; 



2tt?i 2 /? 



where ZNewt. is the one-particle Newtonian partition function. We note that 
a"(0) is a function of x 1 ,x 2 ,x i and it follows from Eq. ((32|) that |ma"(0) may 
be identified as the Newtonian, nonrelativistic potential energy U(x 1 , x 2 , x 3 ) of 
a test particle of mass m with coordinates (0, x 1 , x 2 ,x 3 ) due to the gravitational 
field. This potential energy function U(x 1 , x 2 ,x 3 ) is normalized (by an additive 
constant) so that [7(0,0,0) = 0, i.e., the potential energy on the timelike path 
a is zero. This is because K is the four velocity on a so a(X) = \\K\\/c is iden- 
tically 1 as a function of A on er, forcing a"(0) = there. Examples are given 
in section 7. 

Returning to the general case of an N particle system, let, 

F x (J3,V,N) = -±hiZx, (33) 

where the subscript indicates dependence on the parameter A and Z\ is given 
by Eq. ([M)l . Again, assuming a(0) = 1 and a'(0) = 0, it now follows that, 



liin F\(V,(3,N) = F Ncwt . (V,0,N), (34) 

where, 

F N ew t. 08, V, N) = - ~ In Z Nowt . , (35) 

and 



^Nowt. 



Ill 



M(x- 4 ' w ' ,o, < ix ) ■ w 

The difference between the relativistic and associated Newtonian free energies 
may be computed directly from Eqs. ([24j) and (|36|) . Thus, 



F(f3,V, N) - F Ncwt . (p, V, N) = i In (^^j = jvfcT In (^) 



ni?T In 



/3mc 2 V 2/?mc 2 f ^2l w x 
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where R is the universal gas constant (k times Avogadro's number), and n is 
the number of gram-moles of gas. We note that in the preceding equations, the 
role of volume V is restricted solely to the identification of (constant) limits of 
integration for the volume integrals. In the relativistic context, these limits of 
integration are determined by proper lengths of the dimensions of the container 
of gas, and in the Newtonian limit the limits of integration are absolute length 
measurements . 



7. Equations of State 



In this section we find the equations of state for the relativistic gas particle 
system, and the equation of state for the corresponding Newtonian system of 
particles subject to the potential energy function U(x 1 , x 2 , x 3 ) = ma"(0)/2. For 
notational simplicity and for the purpose of comparison, we use different fonts 
to distinguish the Pressure V and volume V for the Newtonian system from the 
relativistic counterparts P and V. 

Following [9] , the equation of state of the Newtonian particle system, subject 
to the gravitational potential t/(x) = U(x 1 , x 2 , x 3 ) = ma"(0)/2, is given by a 
local version of Boyle's Law, 

PP(X) = PNewt.(x), (38) 

where 'P(x) is local pressure, and /?Ncwt.(x)dx is the probabilistic expectation 
of the number particle^) at position x given by, 

jVexp(-/3t/(x)) 

PNewt.(x = 7 ■ (39) 

J exp(-/3C/(x)))dx 



Integrating both sides of Eq. (|39|) over the volume gives 

(V)V = NkT, (40) 
where (V) = h f v V(x)dx is the averaged pressure. 



An analogous equation of state, 



/3P(x) = p(x), (41) 

holds in the relativistic case. However, in this case, the volume element is 
no longer dx, as in the classical case, but rather dfl = ^/det qijdx, where qij 
is the metric induced by g a on the space slice at fixed time coordinate r = 

3 More formally the integral of PNewt.( x ) over a Borel measurable subset of the region 
containing the gas particles gives the expected number of particles in that set with respect to 
the Canonical Gibbs measure. 
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0, orthogonal to the observer's four-velocity; qtj is the metric that measures 
distances within the container of gas particles. The particle number expectation 
density function, p(x), with respect to the space volume measure dQ is then 
given by, 

„(*) = - gjfefrW) , (42) 

7 (x) v /det 9 y J 1/ (X 2 (7(x))/ 7 (x))dx 



where 7(x) = 7 = a/?mc 2 , as before. Thus, integrating both sides of Eq. (j4"Tjl 
over the volume V with respect to dfl gives, 

{P)V = NkT, (43) 
where (P) — y J v P(x)dfl is the relativistic averaged pressure. 

From Eqs. f4"0)) and (|4"3")l , it follows that at temperature T, the relativistic volume 
and averaged pressure are related to their Newtonian counterparts by, 

V L ldx 

(P) = UV) = J ^__ {V), (44) 

where both integrals are evaluated over the same limits of integration, signified 
by the subscript V as described in the paragraph following Eq. (f37|) . 



As an illustration of Eq. (|44"]) . consider a container of gas in Schwarzschild space- 
time with fixed space coordinates r = r 0} 9 = 7r/2,0 = (see Eq. (foTj) below). 
Using Eq. (53) of [10] . one may readily calculate yj det qtj to third order, and 
with more work to higher order, but for the sake of concision, we display the 
result only to 0(2), 

v/d^7 = 1 + (2 (^) 2 ~ (x 2 ) 2 ~ (x 3 ) 2 ) + 0(3). (45) 

It is apparent that relativistic gravitational effects are negligible unless GM is 
of the order of c 2 or higher. 



8. Examples in Kerr spacetime 



In this section we calculate a(X) and the limiting Newtonian potential energy 
function, U(x 1 1 x 2 , x 3 ) = ma" (0)/2, the key ingredients in Eqs. ([24|) and (|36|) . 
for the cases of circular geodesic orbits in the equatorial plane of Kerr spacetime 
and for a zero angular momentum observer, or ZAMO. 
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The Kerr metric in Boyer-Lindquist coordinates is given by, 

, , / 2X 2 GMr\ 1 l2 4 A 2 G Afar sin 2 6» , £ . 2 

ds 2 = - 1 = — dr = dtd0 + sin 2 

V P J x P P (46) 

+ P-dr 2 + p 2 d6 2 , 
A 

where, 

p 2 = r 2 + A 2 a 2 cos 2 0, 

A = r 2 - 2GM\ 2 r + A 2 a 2 , (47) 
S = (r 2 + A 2 a 2 ) 2 - A 2 a 2 Asin 2 6>, 

and where G is the gravitational constant, M is mass, and a is the angular 
momentum per unit mass, and —GMX ^ a ^ GMX. 

Below we will need notation for A, and E, evaluated at the specific coordinates 
r = ?'o, and 9 = it/2. For that purpose, we define, 

A =rl - 2GMX 2 r + X 2 a 2 , 

9 (48) 

(JZ , \2 2\ 2 \ 2 2 a V ; 

Eo = (r + A a J - A a A . 

Example 1: Circular Geodesic Orbit The circular orbit in the equatorial 
plane (see, e.g., [IT], [12]) with d<j)/dt > at radial coordinate ro, is given by, 



tr(t)= t,r ,7r/2, == , (49) 

When a > 0, this orbit is co-rotational and when a < the orbit is retrograde. 
All circular orbits in the equatorial plane are stable for ro > 9GMA 2 , indepen- 
dent of a, but the co-rotational circular orbit is stable for smaller values, even 
inside the ergosphere at ro = 2GMX 2 when a is sufficiently close to GMX 2 (cf. 

El)- 

The Killing field tangent to the path (where it is the four velocity) is given by, 
/?=( A2 "^ +r ° 3/ % .0,^). (50) 

where, 

V 



\f2X 2 aVGMr ( Q /2) + r% - 3X 2 GMr 2 . (51) 
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Using Eas.lflB]) and ([50]), the calculation for a(A) = X^-K a K a in Boyer- 
Lindquist coordinates, and thence a"(0), is straightforward. The result is, 

1 11, . GMm GMmr 2 sin 2 6» 'iGMm , . 

2 W '(°) = -- 2^ + < 52 > 

This expression may be interpreted as Newtonian potential energy of a single 
particle in a container of particles in circular orbit with angular velocity, 



• # GM 

^^ = V^r (53) 

where ro is the radius of the orbit. The first term on the right hand side of 
Eq. (p)2")) is the gravitational potential energy, and the third term is an additive 
constant that forces a"(0) = when r — ro and 9 = n/2, i.e., at the origin of 
coordinates for the rotating container. The second term is the centrifugal poten- 
tial energy, more readily recognized when expressed in cylindrical coordinates. 
To that end, let p = r sin#, with <p the azimuthal angle, and let z measure linear 
distance along the axis of rotation. The magnitude of angular momentum of 
a uniformly rotating particle of mass m with cylindrical coordinates (p, 0, z) is 
£ = mp 2 (j), and the centrifugal potential energy is then given by, 

£ 2 GMmp 2 GMmr 2 sin 2 6 , . 

(54) 



2mp 2 2r 2 2rg 



The minus signs in Eq. (|54p take into account the direction of force, away from 
the central mass. 

In order to compute a"(0) in Fermi coordinates, we select the following tetrad 
vectors in the tangent space at cr(0). 



eo 



, X 2 aVGM + r i Q 3/2) Q Q VgM \ 



V ' 1 7 v 



ei = ^0,^,0,0 

e 2 = (°>°»^> ) . (55) 
/ A 2 VGM(X 2 a 2 + r 2 - 2X 2 a^/GMr Q ) X 2 VGM(a - 2^) + rf ,2) \ 

e 3 — ^ m— >0,0,- 



This tetrad may be extended via parallel transport to the entire circular orbit 
given by Eq. (|()2"|) , but we need these tetrad vectors only at <r(0). We note that 



1G 



in the case that a = 0, this is the same tetrad utilized in pQ for the estimation 
of the free energy of a gas in circular orbit in Schwarzshild spacetime. 



Fermi coordinates relative to these coordinate axes (i.e., the above tetrad) may 
be calculated using Eq.(27) of [13]. The result for the Boyer-Lindquist coordi- 
nates r and 9 to second order expressed in the Fermi space coordinates x 1 , x 2 ,x 3 
at (x° = r = t = 0) is, 



2 y /7^x 1 x 2 



7T X 

VA^x 1 (r n X 2 GM - X 2 a 2 ) (x 1 ) 

r 1 



r 2r 3 
A (x 2 ) 2 (r„-A 2 GM)(. 3 ) : 



2rg 2r 2 



(56) 



(57) 



Now, calculating a(A) = Av '—K a K a in Boycr-Lindquist coordinates, substitut- 
ing for r and 9 using Eqs. (f56"|) and (|57|) gives, 



GMX 2 (r 2 + 3 a 2 A 2 - 4 aA 2 VnTGM) (x 2 ) 
a(A) - 1 + 2^V 2 



2 



3GMA 2 A (x 1 ) 2 

V 7 - + 0(3) 



(58) 



2r 2 V 2 



Computing the second derivative with respect to A at A = yields, 

GMm GMmx 1 GMm (2 (x 1 ) 2 - (x 2 ) 2 - (x 3 ) 2 ) 



\ma"{Q) = I 



r r 2 2r 3 

3GMm 



' GMm GMmx 1 GMm ((x 1 ) 2 + (x 3 ) 2 ) 



0(3). 



2r ' r 2 2r 3 / ' 2r 

(59) 



Eq.((59|) may be compared term-by-term with Eq. ([52]) . The expression in the 
first pair of parentheses on the right hand side of Eq. (|59|) is the Taylor exapan- 
sion to second order of the gravitational potential, which is the first term on 
the right hand side of Eq. ([52")) . The second terms in both equations are related 
analogously. At the point <r(0) in the orbit, the Cartesian variable x 1 in Eq. ([59| 
measures (Newtonian) distance from the origin of coordinates in the radial di- 
rection away from the central mass, x 2 measures distance in the "z direction" 
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parallel to the axis of rotation, and a; 3 measures distance in the tangential di- 
rection, parallel to the motion of the container of gas in orbit. (However, these 
orientations do not hold at other parts of the orbit since the coordinate axes are 
nonrotating.) Eq. (|59[) may obviously be simplified to yield the potential energy 
function, 

, , , „ 1 „, N ZGmMix 1 ) 2 GmM(x 2 ) 2 

U(x\x 2 ,x 3 ) = ~ma"(0) = -j^J- + -J-J- + 0(3). (60) 

z zr zr 

Combining Eq. lfBT))) with Eqs. (l33|) and (f3"6")) immediately yields the Newtonian 
free energy for a dilute gas in circular orbit around a central mass to which the 
relativistic counterpart Eq. (|33|) (with A = 1/c) may be compared. 

In the special case that a = 0, the Kerr metric (Eq. (|46|) ) reduces to the Schwarzschild 
metric, 

2GM A 2 \ 1 dr 2 

7 ' 



ds 2 = -[l — —dt 2 + 2GMA2N + r 2 (d6 2 + sin 2 Odcfi 1 ). (61) 



In this case, Eq. (f6T)|) may be derived in an alternative manner as follows. The 
Killing field K tangential to the circular geodesic orbit with radial coordinate rg 
was calculated in (T3] to second order in the Fermi coordinates, using the tetrad 
given by Eqs. ([55|) with a = 0. With dependence on G and A made explicit, it 
was found (see Eq.(61) of [13]) that, 



K(r, x\x 2 ,x 3 ) = \l, ^ x\ 0, x 1 J + 0(3). (62) 

Using Eq. (|62p . the Newtonian limit of the energy per particle, 



E = -K a Pa -mc 2 = - (po + ^) + (x 3 Pi - x^) , (63) 

may be computed directly. The limit of the first term on the right side of Eq. (f6"3"]l 
was calculated in [T] as, 

/ mN ( p 2 GMm(2(xi) 2 -(x 2 ) 2 -(x*) 2 ) \ 

A y pa + y) = -[^ 2^ J ' (64) 

where p 2 = p 2 + p\ + p 2 (and because of the use of different unit conventions, 
> " in Eq.jMJ) corresponds to "po/A" in Thus, 
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p 2 GMm (2(,y-(x 2 ) 2 -(x 3 f) /gm 3 ! 
1 ™^ = ^7 53 + \/— Tv* Pi-*Pa)- (65) 



Completing the square for the sum of the first and third terms on the right hand 
side of Eq. (fo5|) gives in the case of p\ , 



p\ t I GM 1 / GM _3l GmM (r 3 ) 

(66) 



2m ' \/ 77S T ^ 1 2m 1 ^ + \/ r„ TOX I 2r^ 



? GmM 



Pi 

2m 2rl 



x 



with a similar expression for the quadratic polynomial in p%. The result is, 



p\ pi p\ ZGmMix 1 ) GmMix 2 ) 

Km E = -p- + p- + , v ; + \ ' + 0(3), 67 

A^o+ 2m 2m 2m 2r% 2rl \ h \ ) 



where, 



GM , 

P3 = P3 - \ — g- mx . (68) 



Since Lebesgue measure dpi is invariant under translations on R, Eq. (f67|) is con- 
sistent with Eq. (f60|) . 

Example 2: Zero Angular Momentum Observer The Killing field for a 
Zero Angular Momentum Observer, or ZAMO, in Boyer-Lindquist coordinates 
in Kerr spacetime is given by, 

- / 1 /XkT „ „ 2GMa\ 2 \ 
Using Eqs.gni) and the calculation for a(X) = X^-K a K a yields, 
' (p 2 - 2GM\ 2 r) E 2 , - 4G 2 A/ 2 a 2 A 6 r (r £ - 2r£ ) sin 2 9 



a ( A ) = 2 2A v 

The limiting Newtonian potential is readily calculated as, 



(70) 
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1 GmM GmM 

-ma"(0) = — — , (71) 

as expected. To find the corresponding expressions in Fermi- Walker coordinates, 
we use the following tetrad of vectors in the tangent space at er(0), 



1 /S n n 2GMa\ 2 
e o = — W-t-,0,0, 



r V Ao ' ' ' j ' 



ei = 



,0,0 



/So 
ro 



= (0,0,-,0). 



e 3 



0,0,0, 



rp 



(72) 

(73) 
(74) 
(75) 



Again using Eq. (27) of [13] we find the transformation from Boyer-Lindquist to 
Fermi- Walker coordinates is given by Eq. (|56p and, 



r = r Q + 



1 

X 



(r GM - a 2 ) A 2 



2r 3 



Ao 

2r 3 



Ao (r 3 



GMa 2 \ 4 



(76) 



2r 2 E 



(, 3 ) 2 + 0(3). 



Now, combining Eqs. (|56")) . ([701) , and ([76]) gives, 



a(x x , x 2 , x 3 , X) = 1 + 



GMX 2 A x GAIX 2 B 1 



■ x 



where 



A D " 2r 4 D 2 
2 B 2 , 2 , 2 GMX 2 j 
2r 4 D {X > + 2r 4 D 2 



(77) 



A = r% + 2r 2 a 2 X 2 + a 4 A 4 - 4GA/r a 2 A 4 , 
B 1 = 2r 7 + 7rla 2 X 2 + 8r 3 a 4 A 4 + llGMr 4 a 2 A 4 + 3r a 6 A 6 
+ UGMr^X 6 + 7GMa 6 X 8 - AG 2 M 2 r a 4 X 8 , 

(78) 

B 2 = r 4 + 4r 2 a 2 A 2 + 3a 4 A 4 - 4GAir a 2 A 4 , 

B 3 = (rg - GMa 2 X 4 ) A, 

D = rl+ r a 2 X 2 + 2GMa 2 X 4 . 
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Calculating a" (A = 0) gives the limiting Newtonian potential energy function, 




GmM 




(x 3 )'+0(3), (79) 



2 



which is the Taylor expansion of the Newtonian potential energy function —GMm/r 
plus an additive constant. 

9. Conclusions 

For spacetimes with timelike Killing fields, we have developed the canonical 
ensemble for an ideal gas whose container is at rest relative to what we call 
"Fermi- Walker-Killing" coordinates. The Hclmholtz free energy for such a sys- 
tem is given by Eqs.([2"i|) and ([25]) . 

The Boltzmann factor in these expressions is invariant with respect to coordinate 
changes, so one may use an arbitrary coordinate system, with the appropriate 
transformation of the volume form for phase space in the integral expressions of 
our formulas. However, the Fermi- Walker-Killing coordinate system is the natu- 
ral choice for statistical mechanics because the container for the particle system 
is at rest in this frame, and these coordinates lead to a general formula for New- 
tonian limits, for the purpose of comparison. Under fairly general conditions 
(namely, a(0) = 1 and a'(0) = 0), the Newtonian limits in these spacetimes are 
given by Eqs. (f55|) and ([5u]) . Because of the inhomogeneity of the geometry of 
space in the relativistic case, and the corresponding non translation invariance 
of the gravitational potential energy function in the Newtonian analogs, pressure 
is necessarily a function of position. The relationship between the relativistic 
and Newtonian pressures, averaged over volumes, is given by Eq. (|44p . 

The generalization of our results to a gas with interactions might be carried out 
by starting with a covariant Lagrangian that incorporates the interaction be- 
tween the particles in the container. An energy momentum tensor determined 
by such a Lagrangian together with the Einstein field equations then in principle 
determines the metric. From there the metric together with a timelike Killing 
field (or perhaps an "approximate Killing field"), following the formalism of this 
paper, may be used to develop the canonical ensemble. A simpler, approximate 
theory for nonrelativistic particles with interactions, for the purpose of studying 
the effects of the background curvature could be achieved, following the classi- 
cal formulations, by including a potential energy function of the proper distance 
between particles in the Boltzmann factor in the partition function. 

However, an important physical observation is already available from the re- 
sults of this paper for noninteracting particles. As noted earlier, relativistic 
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corrections to the pressure given by Eq. (|4"4"|) become significant (at arbitrary 
temperature) only in the presence of extremely large masses. It is well-known 
that special relativistic corrections to the Helmholtz free energy for the ideal 
gas are insignificant except at extremely high temperature [7J. It follows from 
our results for the general relativistic case, in particular Eq. ([3"T|) . that curvature 
corrections from strong gravitational fields (so that a varies rapidly as a func- 
tion of position) to the free energy become significant only at temperatures of 
the same magnitude as in the special relativistic case. 
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